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Abstract 

Low order equivalent system models were 
identified from flight test data for the Tu-144LL 
supersonic transport aircraft. Flight test maneuvers 
were executed by Russian and American test pilots 
Hying the aircraft from Zhukovsky airfield outside 
Moscow, Russia. Flight tests included longitudinal and 
lateral / directional maneuvers at supersonic cruise 
flight conditions. Piloted frequency sweeps and 
multi-step maneuvers were used to generate data for 
closed loop low order equivalent system modeling. 
Model parameters were estimated using a flexible, high 
accuracy Fourier transform and an equation error/ 
output error (EE/OE) formulation in the frequency 
domain. Results were compared to parameter 
estimates obtained using spectral estimation and 
subsequent least squares fit to frequency response data 
in Bode plots. Modeling results from the two methods 
agreed well for both a frequency sweep and multiple 
concatenated multi-step maneuvers. For a single 
multi-step maneuver, the EE/OE method gave a better 
model fit with improved prediction capability. A 
summary of closed loop low order equivalent system 
identification results for the Tu-144LL, including 
estimated parameters, standard errors, and flying 
qualities level predictions, were computed and 
tabulated. 

Nomenclature 

e equation error 

EXT extended 

FQ Flying Qualities 
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h altitude, ft 

HSR High Speed Research 

I x ,l Y ,f : body-axis moments of inertia, slug-ft“ 
j imaginary number V— T 

J cost function 

LOES Low Order Equivalent System 

M Mach number 

rn number of discrete frequencies 

p. c/, r roll, pitch, and yaw rates, deg/sec 
RET retracted 

.v Laplace transform variable 

t time, sec 

Tft rol) mode time constant, sec 

u input 

v output error 

w weight, lbf 

_V ( .^ longitudinal center of gravity position 

y model output 

z measured output 

a angle of attack, deg 

tj pilot control deflection, cm 

r equivalent time delay, sec 

co frequency, rad/sec 

tf damping ratio 

Z phase angle, deg 

| | magnitude of a complex number 

superscripts 

f complex conjugate transpose 

estimate 

Laplace or Fourier transform 
-1 matrix inverse 
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subscripts 


a 

lateral stick 

DR 

Dutch Roll 

e 

longitudinal stick 

EE 

Equation Error 

LOES 

Low Order Equivalent System 

LS 

Least Squares 

OE 

Output Error 

0 

nominal or trim value 

r 

rudder pedal 

SP 

Short Period 


Introduction 

As part of the NASA High Speed Research 
(HSR) program, the Tupolev Aircraft Company in 
Russia was commissioned to refurbish, re-engine, and 
instrument an out-of-service Tu-144D supersonic 
transport aircraft, for use as a supersonic flying 
laboratory. The modified aircraft, designated the 
Tu-144LL, was operated from Zhukovsky airfield, near 
Moscow, Russia. 

Flight test data was collected for 6 experimental 
investigations related to high speed research. One of 
the experimental investigations was closed loop flying 
qualities. Flight testing began with 19 Phase I flights 
in 1996-1997, flown exclusively by Tupolev pilots. 

Ref. [1] contains results of longitudinal, lateral, 
and directional closed loop Low Order Equivalent 
System (LOES) modeling based on the Phase 1 flight 
test data, using an equation error method in the 
frequency domain 1 \ The LOES model characterizes 
the linearized dynamic closed loop response of the 
airframe and control system as it appears to the pilot. 
The form of the LOES model is the same as that for an 
open loop unaugmented airplane with classical 
dynamic modes, except the inputs are pilot controls 
with equivalent time delays, instead of control surface 
deflections. Equivalent time delay is introduced to 
account for time delay resulting from the control 
system implementation (e.g., sampling delay), and the 
phase lag at high frequency from control system 
dynamics and various nonlinearities, such as control 
surface rate limiting. 

In the past, the majority of flight testing for flying 
qualities evaluation and specification compliance was 
done using unaugmented aircraft with classical 
dynamic modal responses. To apply specifications 
developed for unaugmented airplanes to airplanes 


whose dynamics are no longer classical because of 
additional dynamics from the control system, the 
concept of a Low Order Equivalent System (LOES) 
model was introduced 4-6 . 

If the closed loop response of an augmented 
airplane to pilot inputs can be accurately characterized 
by a LOES model, then specifications for classical 
unaugmented aircraft parameters found in the current 
Military Specification for Flying Qualities of Piloted 
Airplanes 4 (hereafter called the Mil-Spec) can be 
applied directly to the estimated parameters in the 
LOES model. The Mil-Spec quantifies the relationship 
between parameter values in a low order dynamic 
model and pilot opinion, as measured by flying 
qualities levels 4 , which are based on Cooper-Harper 
ratings 6 . The Mil-Spec also contains specifications 
relating flying qualities levels and equivalent time 
delay. LOES models for the aircraft closed loop 
dynamic response are fit over the frequency band 
corresponding to typical pilot inputs, 0. 1 - 10 rad/sec. 

Many flight test research programs 4 - 6 " 9 have 
demonstrated that the LOES concept can be used to 
predict pilot opinion (quantified by flying qualities 
level) for augmented aircraft dynamic response that is 
in reality high order and/or nonlinear. LOES models 
identified from flight test data are also useful for 
validating linear control law designs and verifying 
specification compliance, since the LOES model 
represents the achieved linearized closed loop 
dynamics of the aircraft. LOES models can also be 
used for rudimentary simulation within the limited 
flight envelope where the model is valid. 

A second phase of flight testing for closed loop 
flying qualities research was conducted on the 
Tu-144LL in the fall of 1998. Phase II testing included 
four flights - one subsonic flight with an all Russian 
crew (flight 20), one subsonic flight with an American 
evaluation pilot (flight 21), and two supersonic flights 
with American evaluation pilots (flights 22 and 23). 
These flights included flying qualities evaluation tasks, 
approaches, and special maneuvers designed to collect 
data for closed loop LOES modeling. Ref. [10] 
contains extensive information on the Tu-144LL 
aircraft, flight profiles, and maneuvers for these flights, 
including qualitative pilot commentary and narrative 
flight summaries written by the American evaluation 
pilots. 

This work examined LOES modeling based on 
the Tu-144LL Phase II flight test data using an 
equation error / output error technique in the frequency 
domain with a flexible, high accuracy Fourier 
transform 2 * . Results were compared with LOES models 
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obtained from a least squares fit to an empirical 
estimate of the Bode plot, using commercially 
available software called CIFER® 1 Correlation of 
the analytical LOES modeling results with pilot 
opinion was not done, because the pilots gave only 
qualitative narrative assessments, and not 
Cooper-Harper ratings. 

The next section describes the two LOES model 
parameter estimation techniques used. Next, LOES 
modeling issues are examined by applying the two 
parameter estimation techniques to Tu-144LL Phase II 
flight test data. Finally, LOES modeling results for the 
Tu-144LL Phase II flights are computed and tabulated. 

Methods 

The model structure for LOES modeling is fixed 
a priori to correspond to classical linear aircraft 
dynamics with an input time delay. For the short 
period longitudinal dynamic mode, the closed loop 
pitch rate response to longitudinal stick deflection is 
modeled in transfer function form as 4 


The forms of the models in Eqs. (2) and (3) follow 
from the closed loop lateral / directional state space 
model, omitting the gravity term in the side force 
equation. 

The problem addressed in this work is accurate 
estimation of the model parameters in Eqs. ( 1 )-(3) 
using measured input-output flight test data. The idea 
is to match the measured outputs or output time 
derivatives with the corresponding quantities from the 
model by adjusting model parameters to minimize a 
measure of fit error, which is the sum of the squared 
deviations between model quantity and measured 
quantity. Several methods exist for estimating model 
parameters based on measured data, both in the time 
domain 14 ’ 15 and the frequency domain 1 ' 4 * 1 M \ 

For identification purposes, Eqs. (1 M3) can be 
re-parameterized as 

- {l’\s+b 0 )e~ T ‘- x 

q = — — T) e (4) 
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The current Mil-Spec 4 correlates pilot opinion, 
quantified by flying qualities levels from the 
Cooper-Harper scale 5 , with ranges of values for all 
model parameters in Eq. (1). except Kq . 

For lateral / directional closed loop dynamics, the 
LOES model for roll rate response is 
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For closed loop yaw rate response, the LOES model is 
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The relationship among the parameters in 
Eqs. (I )-(3) and Eqs. (4)-(6) can be found by direct 
comparison. All LOES modeling was done in the 
frequency domain, because the equivalent time delay 
parameters can be estimated more easily and accurately 
in the frequency domain 24 *. 

Two methods were used to estimate the 
parameters in Eqs. (4)-(6) based on flight test data. 

The first method, described in Refs. [ 1 1 ]-[ 13], involves 
computing a non-parametric estimate of the frequency 
response in the form of a Bode plot. Then, the 
parameters in the models of Eqs. (4 )-( 6) are adjusted to 
minimize the least squares cost function: 
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Typically, the number of frequencies m is 50, and 
the frequencies are chosen logarithmically on the 
interval [0.1,10] rad/sec. Commercially available 
FORTRAN software called CIFER® 1 113 was used to 
estimate LOES model parameters with this method. 

The second method, described in Ref. [2], is a 
two-step approach using equation error (EE) and output 
error (OE) formulations in the frequency domain for 
the parameter estimation. For the equation error 
formulation, the cost function is: 


I '» | 

J EE = -£ D (M) ~ MM)«/ € ~ JC ° iX 
Z /=l 




(8) 


where N(jco } ) and D(jco i ) are the numerator and 
denominator polynomials of the transfer function, with 
s = jo)f . Accurate Fourier transforms for m arbitrary 
frequencies were computed using the method of 
Ref. [16]. Parameter estimation results from the 
equation error method were used as starting values for 
the output error method, which minimized the cost 
function: 


t m i w . i 

Joe ! =^W) (91 

Final parameter estimates were from the output 
error formulation, because of known favorable 
asymptotic properties 14 ' 15 . Equation error is much 
more robust to starting values of the parameters, 
because the model output depends linearly on all model 
parameters except the equivalent time delay T. In this 
work, the equation error method always converged 
quickly and without difficulty using zeros for the 
starting values of all the parameters. Parameter 
estimates from the equation error method provide good 
starting values for the output error method. The output 
error method needs good starting values for the 
parameters because the model output depends 
nonlinearly on the transfer function denominator 
parameters, in addition to the equivalent time delay r. 


The EE/OE sequence therefore gives parameter 
estimation results with good asymptotic properties, 
without requiring good starting values for the 
parameters. Estimated parameter standard errors do 
not require correction, because the analysis in the 
frequency domain automatically accounts for colored 
residuals 2 ’ 14,15 . 

Both the equation error and output error methods 
are nonlinear estimation problems, because of the 
equivalent time delay parameter. Details on the 
optimization methods used to find both equation error 
and output error parameter estimates can be found in 
Ref. [2]. 

Any of the methods described above can be used 
for general multi-input, multi-output models. For the 
results given here, Eqs. (4)-(6) were used to identify 
LOES models from flight test data. 

Important differences in the two methods can be 
appreciated by examining the cost functions in 
Eqs. (7)-(9). The CIFER® method minimizes a cost 
comprised of two parts - the squared difference 
between measured and model transfer function 
frequency response magnitude in dB and phase in deg. 
The first step in the CIFER® method is non-parametric 
estimation of the transfer function frequency response 
in the form of a Bode plot, obtained using the ratio of 
estimates for the input-output cross-spectral density 
and input auto-spectral density. The two weighted 
terms in the cost involve nonlinear functions of the 
complex frequency response (gain in dB and phase 
angle in deg). Since the frequency response estimate is 
calculated from a ratio of spectral estimates, the 
accuracy of the estimated frequency response can 
degrade when the frequency content of the data is 
sparse. CIFER® uses coherence weighting in the 
parameter estimation to emphasize frequency domain 
data with good linearity and signal-to-noise ratio. The 
EE/OE method uses high accuracy, arbitrary frequency 
Fourier transforms 16 directly, and minimizes the square 
magnitude of the equation error (EE) or the output 
error (OE), with equal weighting of the real and 
imaginary parts of the complex quantities in the cost. 
Spectral components in the data, whether sparse or not, 
influence parameter estimation results according to 
their complex magnitude in the frequency domain. 

Details of the calculations and optimization 
methods used to find the parameter estimates that 
minimized the cost functions in Eqs. (7)-(9) can be 
found in Refs. [ 1 1 ]-[ 1 3] for the CIFER® method, and 
Ref. [2] for the EE/OE method. 
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Aircraft 

The Tu-144LL test aircraft is a modified version 
of the Tupolev Tu-144D cranked delta wing supersonic 
transport. Modifications included installation of 
Kuznetsov NK-32 1 engines rated at 3 1 ,000 Ibf static 
dry thrust and 55,000 lbf in afterburner at sea level, and 
installation of flight-test quality instrumentation with 
sampling rate 32 Hz. A 3-view drawing of the 
Tu-144LL appears as Figure 1. Vehicle characteristics 
are summarized in Table 1 . For cruise flight, the 
canards are retracted along the fuselage above and 
behind the cockpit, and the nose is raised. On 
approach, the canards are extended for slow flight and 
the nose is drooped for enhanced visibility from the 
cockpit. Flight crew consisted of two pilots, a 
navigator, and a flight engineer. Details of the flight 
control system are available in Ref. [ 10|. 

Air data on the Tu- 144LL was obtained from 
production sensors, not specialized flight test 
instrumentation, such as an air data probe. Data 
compatibility analysis 17 indicated that measurements 
of angle of attack and sideslip angle contained large 
systematic errors, including lag and hysteresis effects. 
Consequently, the flight test data analysis was done 
using only body-axis angular rate data for the measured 
outputs. The angular rate data was found to be 
consistent with measured Euler angle data, using data 
compatibility analysis with the rotational kinematic 
equations only 17 . Finally, there was some evidence to 
suggest that time skews were present in the measured 
input-output data, due to the instrumentation system. 
The main clue was an excessive time delay 
(-200 msec) between initial movement of control 
surfaces and angular rate response. However, it could 
not be determined conclusively what, if anything, the 
time skews were, so the data was analyzed as recorded. 

Results 

The first case studied was closed loop 
identification of the Tu-144LL longitudinal short 
period LOES dynamics for supersonic cruise flight at 
Mach 1.93. The pilot applied a frequency sweep input 
to the longitudinal stick, shown in the upper plot of 
Figure 2. Measured pitch rate response is shown as the 
solid line in the middle plot of Figure 2. The dotted 
line on the same plot shows the computed pitch rate 
response to measured stick input, using the LOES 
model of Eq. (4) with parameters estimated using the 
EE/OE method. The lower plot in Figure 2 shows the 
output residual, which is the difference between the 
measured and computed pitch rate traces shown in the 


middle plot of Figure 2. Frequencies used in the data 
analysis were evenly spaced at 0. 1 rad/sec on the 
interval [0.1,10] rad/sec. The model fit in the 
frequency domain (not shown) was excellent, and 
Figure 2 shows that the model fit in the time domain 
was also excellent. Estimated parameter values and 
standard errors for the identified q /q e LOES model 
are listed in column 2 of Table 2. Estimated standard 
errors in Table 2 indicate that the model parameter 
accuracy was better than 8%, based on parameter 
estimate values, for all q jl] e LOES model parameters. 

CIFER® was used with the same data to estimate 
model parameters in the same LOES model of Eq. (4). 
Results are given in the third column of Table 2. 
Comparison of the values in columns 2 and 3 of 
Table 2 indicates that the results from CIFER® agreed 
well with corresponding values from the EE/OE 
method. This built confidence in the results, because 
the two methods differ in formulation, cost function, 
and optimization technique, but gave similar LOES 
modeling results. Since the CIFER® model parameters 
were so close to the EE/OE model parameters, the 
computed pitch rate response from the CIFER® model 
was nearly the same as the computed pitch rate 
response from the EE/OE model, shown by the dotted 
line in Figure 2. 

The upper plot of Figure 3 shows a different kind 
of longitudinal stick input implemented by the pilot at 
approximately the same flight condition. This input 
was a modified doublet called a 2- 1 - 1 , which ideally 
consists of alternating pulses with pulse widths in the 
ratio 2- 1 - 1 . Inputs of this type are also called 
multi-step inputs. The widths of the pulses were 
chosen so that the frequencies corresponding to the “2” 
and “1” pulses bracketed the a priori estimate of the 
natural frequency for the closed loop short period 
mode. Data analysis results from Phase I nights 1 at 
transonic cruise conditions were used to design the 
input pulse width and amplitude. 

The EE/OE and CIFER® methods were used to 
estimate LOES model parameters in Eq. (4), based on 
data from the 2-1-1 maneuver. Table 3 contains the 
results. The EE/OE method produced LOES model 
parameter estimates close to those estimated using the 
frequency sweep input at approximately the same flight 
condition (cf. the second column of Tables 2 and 3). 
Some discrepancies in the parameter estimates may be 
due to slight differences in flight condition and mass 
properties. The lower plot in Figure 3 shows a good 
time domain match to the measured pitch rale for the 
EE/OE model. CIFER® parameter estimates in column 
3 of Table 3 were different than the EE/OE values. 
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mainly because the 2-1-1 data did not provide 
sufficient power across the LOES modeling frequency 
range to allow an accurate Bode plot estimate from 
spectral estimates. Figure 3 shows the corresponding 
time domain mismatch to the pitch rate measurement 
for the CIFER® model based on data from 2-1-1 
maneuver 22_4f only, denoted by CIFER 4f. 

The CIFER® method was not designed to work 
with data of this type. However, results from CIFER® 
improved and approached the EE/OE results when 
additional 2-1-1 maneuvers at the same flight condition 
were concatenated for the analysis, allowing sufficient 
averaging from the data windows for accurate spectral 
estimation. Additional 2-1 -I maneuvers 22_4a and 
22_4b were concatenated to the original single 
maneuver 22_4f for another CIFER® analysis. The 
results shown in column 4 of Table 3 are in good 
agreement with the EE/OE results, and Figure 3 shows 
a good match to the pitch rate measurement. The 
multiple-maneuver CIFER® results are denoted by 
CIFER 4abf. The CIFER® method therefore works 
well for 2-1-1 maneuvers also, given sufficient data for 
good spectral estimates. The EE/OE method is less 
sensitive to the spectral content in the data, because 
Fourier transforms are used directly, instead of spectral 
estimates. This allowed the EE/OE method to 
accurately identify LOES models using data from 
either a frequency sweep or a single 2-1-1 maneuver. 

The EE/OE method can also be used with 
multiple maneuvers. For this case, results from the 
EE/OE method using concatenated data from 
maneuvers 22_4a, 22_4b, and 22_4f were similar to 
those given in column 2 of Table 3, which were based 
on data from maneuver 22_4f only. 

Flight time required for the 2-1-1 maneuver is 
approximately 1/6 the time required for the frequency 
sweep maneuver. It might be possible to shorten the 
frequency sweep maneuver and maintain its 
effectiveness by omitting the lower frequency portions 
of the maneuver, since only the relatively high 
frequency closed loop short period dynamics were 
being identified. 

Figure 4 shows a similar 2-1-1 maneuver 
executed at approximately the same flight condition. 

In the lower plot of Figure 4, the LOES models 
identified in the previous example (cf. Table 3) were 
used with the input shown in the top of Figure 4 to 
predict the pitch rate response. For this prediction 
case, the fidelity of the EE/OE and CIFER® models 
was similar to what was seen in Figure 3. 


Figure 5 shows another longitudinal 2-1-1 
maneuver, this time at transonic cruise, Mach 0.88. 

The model trace shown in the lower plot of Figure 5 
was generated using the LOES model of Eq. (4) with 
parameters estimated using the EE/OE method. 

Results for this case are given in column 3 of Table 4. 

In Figure 6, the identified LOES model was used for a 
prediction case at the same flight condition, with 
inverted polarity of the 2-1-1 longitudinal stick input. 
Figure 6 demonstrates the excellent prediction 
capability for LOES models identified using the EE/OE 
method, and also indicates that the LOES model was 
not sensitive to input polarity for this flight condition. 

The final example was multiple-input, multiple- 
output lateral / directional LOES modeling at 
supersonic cruise. Roll rate and yaw rate outputs were 
used with both lateral stick and rudder pedal inputs, in 
the model structure of Eqs. (5) and (6). Data from both 
a rudder pedal sweep and a lateral stick sweep 
maneuver were combined in the frequency domain for 
LOES modeling using the EE/OE method. The input- 
output data in Figure 7 is from the lateral stick sweep 
maneuver. The model fit to the data in the time 
domain was excellent, demonstrating that the EE/OE 
method can be successfully used for multiple-input, 
multiple-output lateral / directional LOES modeling. 
The quality of the time domain fit to the rudder pedal 
sweep maneuver (not shown) was similar. The last 
column of Table 4 contains results for this case. 

Longitudinal and lateral / directional LOES 
modeling results for Tu-144LL Phase II flights are 
summarized in Table 4. Flight conditions included are: 
canard-extended approach at Mach 0.34, transonic 
cruise at Mach 0.88, and supersonic cruise at Mach 
1 .93. All LOES modeling results in Table 4 were 
obtained using the EE/OE method. More than one 
maneuver was used in cases where the flight conditions 
for the maneuvers were nearly the same. This is 
indicated by the maneuver(s) listed for each flight 
condition. Flying qualities level predictions using the 
Mil-Spec were level 2 or 3 because of large values for 
estimated equivalent time delay, cf. Table 4. Based on 
all LOES parameter estimates except the equivalent 
time delay, the Mil-Spec indicated that the aircraft had 
level 2 longitudinal flying qualities for supersonic 
cruise, transonic cruise, and canard-extended approach, 
and level 1 lateral / directional flying qualities for 
supersonic cruise. Flying qualities levels predicted 
exclusive of the equivalent time delay (i.e., level 2 
longitudinal, and level 1 lateral / directional) were 
consistent with qualitative comments from the 
evaluation pilots 10 . This was further evidence of the 
possible instrumentation time skew mentioned earlier. 
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Similar findings of excessive time delay were reported 
in Ref [ 1 ] for Phase I flight test data analysis. 

LOES modeling was also done using the EE/OE 
method for several longitudinal cases with 150 msec 
delay removed from the pitch rate measurement. This 
was the estimated time delay that could he reasonably 
attributed to instrumentation time skew r based on 
control surface deflection and pilch rate measurements 
for the multi-step maneuvers, allowing for normal pitch 
rate response delay. For pitch rate data with 150 msec 
time delay removed, estimated LOES model 
parameters were essentially unchanged, except that the 
equivalent time delay estimate was roughly 150 msec 
smaller. If estimated equivalent time delays were 
reduced by 150 msec for all longitudinal and lateral / 
directional cases in Table 4, the Mil-Spec indicated that 
the aircraft would have level 2 longitudinal flying 
qualities and level 1 lateral / directional flying 
qualities, again in concurrence with qualitative pilot 
opinion 10 . 

Concluding Remarks 

This work was concerned with accurately 
identifying closed loop Low Order Equivalent System 
(LOES) models from flight test data for the Tu-144LL 
supersonic transport aircraft. LOES models are linear 
models with input time delay, used to characterize 
closed loop aircraft response to pilot inputs. Data 
analysis was done in the frequency domain using two 
methods: a two-step equation error / output error 
technique (EE/OE) using a high accuracy Fourier 
transform with selectable frequency range and 
resolution, and a least squares fit to an empirical 
estimate of the Bode plot (CIFER® ). Both methods 
identified accurate LOES models with good prediction 
capability when used with data from a frequency sweep 
maneuver or multiple concatenated multi-step 
maneuvers. For a single multi-step maneuver, the 
EE/OE method gave a better model fit with improved 
prediction capability. 

The EE/OE modeling technique w-as applied to 
Phase II flight test data from the Tu-144LL aircraft to 
identify closed loop LOES models at various flight 
conditions, including transonic cruise, canard-extended 
approach, and supersonic cruise. Longitudinal and 
lateral / directional LOES modeling results were 
computed and presented. Identified LOES models 
were consistent with qualitative pilot opinion, if 
equivalent time delays were adjusted for estimated 
instrumentation time skew. 


The results shown here suggest that the EE/OE 
method can be used in conjunction with multi-step 
maneuvers to decrease flight test time required to 
collect data for LOES model identification. All LOES 
modeling was done using only measurements for pilot 
inputs and body-axis angular rate outputs, 
demonstrating that this type of analysis can be done 
successfully with minimal instrumentation. 
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Table 1 Tu- I44LL Vehicle Characteristics 


Length, ft 

215.5 

Height, wheels up, ft 

42.2 

Wing Span, ft 

94.5 

Wing Area, ft 2 

6769 

Mean Aerodynamic Chord, ft 

76.5 

Max. Fuel, Ibm 

209,440 

Est. Max. Range, nm 

1,620 

Max. Mach Number 

2.4 


Table 2 Tu-144LL Longitudinal Short Period LOES 
Modeling Results, Frequency Sweep Maneuver 23_la, 


a Q - 4.8 deg, M 0 = 1.93, x CH = 0.46c, 7 V = 1.44xl0 7 


Maneuver 23_ 1 a 

Parameter 

EE/OE 

CIFER® 


Estimate 

Estimate 


(Std. Error) 

(Std. Error) 

ht 

0.353 

0.380 

u \ 

(0.011) 

(0.022) 

*0 

0.106 

0.107 

(0.006) 

(0.013) 

«! 

0.932 

1.040 

(0.035) 

(0.121) 


1.970 

2.103 

a 0 

(0.035) 

(0.133) 

T 

0.194 

0.192 

1 e 

(0.014) 

(0.013) 

'/Te 2 

0.30 

0.28 

Csp 

0.33 

0.36 

co SP 

1.40 

1.45 


Table 3 Tu-144LL Longitudinal Short Period LOES 
Modeling Results, 2-1-1 Maneuvers 


a 0 = 5 . 1 deg, M 0 = 1.90, x Cf , = 0.46c , / v = 1.41 xlO 7 



Maneuver 22_4f 

Maneuvers 

22_4a,b,f 

Parameter 

EE/OE 

CIFER® 

CIFER® 


Estimate 
(Std. Error) 

Estimate 

Estimate 

b \ 

0.375 

(0.014) 

0.304 

0.430 

b o 

0.1 14 
(0.008) 

0.506 

0.170 

“l 

0.905 

(0.032) 

1.046 

1.064 

a 0 

1.746 

(0.031) 

3.337 

1.760 

r. 

0.188 

(0.024) 

0.138 

0.174 

'/r* 

0.30 

1.67 

0.40 

Csp 

0.34 

0.29 

0.40 

m SP 

1.32 

1.83 

1.33 
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Table 4 Tu-144LL LOES Modeling Summary Results - Phase II Flight Tests 



Longitudinal ]j 

| Lateral / Directional 


Supersonic 

Cruise 

Transonic 

Cruise 

Canard 

Extended 

Approach 


Supersonic 

Cruise 

Maneuver(s) 

f23_ 1 a 

f20_4b,c,d 

f20_4h.i,j 

Maneuver(s) 

f23_2a.3a 

Mo 

1.93 

0.88 

0.34 

M„ 

1.94 

a„ (deg) 

4.8 

5.8 

8.1 

a„ (deg) 

4.6 

h„ (KP ft) 

53.4 

29.7 

7.0 

h ( , (10 3 ft) 

53.9 

w (Ibf) 

3.35 x 10 s 

3.22 x I0- S 

3.04 x 10 5 

w (Ibf) 

3.13 x 10 s 

/y (SlUg-ft 2 ) 

1.44 x 1() 7 

1.37 x 10 7 

1.28 x 10 7 

/y (SlUg-ft 2 ) 

1.17 x 10 6 


0.46 c 

0.46 c 

0.41 c 

/, (slug-fl 2 ) 

1.51 x 10 7 

Landing Gear 

RET 

RET 

EXT 

l xz (slug-ft 2 ) 

-9.75 x 1 0 4 

Canard 

RET 

RET 

EXT 

*CK 

0.46 c 

Nose (deg) 

0 

0 

17 

Landing Gear 

RET 


0.353 ±0.011 

1.512 ±0.080 

0.626 ± 0.026 

Canard 

RET 

I’o 

0.1 06 ±0.006 

1.894 ±0.101 

0.457 ±0.1 54 

Nose (deg) 

0 

a\ 

0.932 ± 0.035 

3.680 ±0.1 82 

1.741 ±0.156 

Ci 


a {) 

1 .970 ± 0.035 

7.246 ± 0.286 

1.518 ±0.380 

c \ 



0.194 ±0.014 

0.280 ± 0.009 

0.200 ± 0.008 

c () 


«6 

0.35 

1.51 

0.63 

di 


1/^2 

0.30 

1.25 

0.73 

d \ 

0.490 ±0.1 84 

Csp 

0.33 

0.68 

0.71 

d o 

2.377 ±0.158 

OJ S p 

1.40 

2.69 

1.23 

e 2 

9.514 ±0.521 





n 

1.052 ±0.1 13 





e o 

12.676 ±0.7 17 





f 2 

0.278 ± 0.060 





f i 

1.107 ±0.068 





/() 

0.263 ± 0.033 





82 






#1 

nag 





£o 

1 .082 ± 0.084 






0.221 ±0.039 






0.218 ±0.010 





Cdh 

0.169 





10 DR 

1.27 





l/7> 

4.34 

Mil-Std 1797 A 
FQ Level 

2 

3 

2 

Mil-Std 1797 A 
FQ Level 

3 
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Tu 1 44LL 3- vi ew 

Figure 1 Tu-144LL 3-view 
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Figure 2 Frequency Sweep Maneuver 23_ 1 a, 
a a = 4.8 deg, M () = 1 .93, x CK = 0.46c , / v = 1 .44x 1 0 7 


flight 



Figure 3 2-1-1 Maneuver 22_4f, 
a G - 5. 1 deg, M () = 1 .90, x cg = 0.46F , / v = 1 .4 1 x 1 0 7 



Figure 4 2-1-1 Prediction Maneuver 22_4a, 
a G = 5.1 deg, M a = 1 .9 1 , jr r# = 0.46c , / v = 1.44x1 0 7 


10 

American Institute of Aeronautics and Astronautics 







Figure 5 2-1-1 Identification Maneuver 20_4d, 

= 5.8 deg, M a = 0.88, x cx - 0.46c , / v = 1.37x1 0 7 




Figure 6 2-1-1 Prediction Maneuver 20_4f, 

= 5.8 deg, M () = 0.88, x CM = 0.46c, / v = 1.37xl0 7 
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Figure 7 Frequency Sweep Maneuver 23_2a, 
a () = 4.6 deg, M t> = 1 .94, x cg = 0.46c . 

/, = 1.17x1 0 6 , /. = 1.51x1 0 7 . I x: =-9.75xl0 4 
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